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Abstract 

We calculate the mixed tensor susceptibility of QCD vacuum in the framework of the 
global color symmetry model. In our calculation, the functional integration over gluon fields 
can be performed and the gluonic vacuum observable can be expressed in terms of the quark 
operators and the gluon propagator. The mixed tensor susceptibility was obtained with the 
subtraction of the perturbative contribution which is evaluated by the Wigner solution of 
the quark gap equation. Using several different effective quark-quark interaction models, 
we find the values of the mixed tensor susceptibility are very small. 
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1. Introduction 

In SVZ sum rules, in order to determine the static properties of hadrons it was suggested 
to consider two-point correlator functions of quark currents in the presence of an external 
constant classical field, where nonperturbative effects are taken into account in the so-called 
vacuum susceptibilities pQ. These induced condensates play important roles in determination 
of the hadron properties such as the nucleon magnetic moments^, the isovector axial coupling 
constant El E] j the isoscalar axial coupling constant 4., the pion- nucleon coupling constant [H] 
and the nucleon tensor charge 6, Zj within this version of SVZ sum rules. In the literature, 
there are always two kinds of vacuum susceptibility that appear in the conventional two-point 
treatment of an external current field : one is the induced quark condensate and the other 
is the induced mixed quark-gluon condensate. For convenience, we refer to the former as the 
quark condensate susceptibility and the later as mixed condensate susceptibility in this letter. 

The vacuum tensor susceptibilities are relevant for the determination of nucleon tensor 
charge^ E|. The value of nucleon charge is related to the first moment of the transversity 
distribution h±(x), where h±(x) is an additional twist-two chirality violating structure function 
which can be measured in the Drell-Yan process with both beam and target transversely polar- 
ized. The previous evaluation of the quark condensate tensor susceptibility were performed in 
the framework of QCD sum rules |S1 El El El- the chiral constituent model^U] and global color 
symmetry model(GCM) [XT] respectively. Actually, there still exist uncertainty about this in- 
duced susceptibility since different theoretical treatments can give very different results, which 
should be checked by the future measurement of the transversity distribution h\(x). Another 
tensor susceptibility, the mixed condensate tensor susceptibility was only evaluated roughly 
in Ref 6 . within the two-point function of QCD sum rules. In this letter, we will give the 
calculation of the vacuum mixed tensor susceptibility within the framework of GCM in the 
mean field level. 

As a truncated DSE-model, GCM is a quite successful four-fermion interaction field theory 
which can be directly derived through a truncation of QCD l21E]- This truncated DSE model 
has been applied extensively at zero temperature and chemical potential to the phenomenology 
of QCD [T5lll6j . including the studies of observables from strong interaction to weak interaction 
area. Furthermore, the truncated DSE models also have made important progress in the studies 
of strong QCD at finite temperature and chemical potential JJj. Due to the fact that the 
evaluation of mixed quark-gluon condensate can be performed in the framework of GCM[T2*]. 
a method to evaluate the mixed vacuum tensor susceptibility is proposed within this DSE 
formalism and the numerical results for the mixed tensor susceptibility are given. 

2. Formalism 

In the chiral limit, the QCD generating functional for quark field in the Euclidean space is 
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given by 



where 



Z[r], rj\ = j VqVqVw'Dw'DAexp / — S — S g f — S g + J d 4 x(rjq + 577) | , 
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and 5 9 /, S g are the gauge-fixing and ghost actions respectively. Through introducing a bilocal 
field B e {x,y) as in the generating functional can be given as 
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with 



Sfo^B^y)] = j j d A xd%\q{x)G~ l {x,y-[B d \)q{y) 
, B B {x,y)B 6 {y,x)' 



2g 2 D(x - y) 



(5) 



G-V, y; = 7 • - y) + A 9 5 e (x, y). (6) 

The quantity A 61 arises from Fierz reordering transformation and is the direct product of Dirac, 
flavor SU(3) and color matrices 



(7) 



And g 2 D(x — y) is the connected gluon two-point function without quark-loop contributions. 
Neglecting W\ [J^], we can obtain the GCM generating functional as 



Zgcm[v,v] = J T ) q'DqDB (x,y)exp^-S[q,q,B (x,y)}+ J d 4 x(i]q + qrj) 



(8) 



Performing the functional integration over T>q and T>q in above equation, we obtain 

ZgcaMv] = J VB e (x,y)eM-S[f],V,B e (x 1 y)}), (9) 
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where 



S[rj,rj,B e (x,y)] = -Trhx d ■ -y6(x - y) + A e B (x,y) 
+ 



B 6 (x,y)B e (y,x) e 

+ r]{x)G{x,y;B )r?(y) 



2g 2 D(x - y) 



(10) 



The saddle point of this action is defined as 8S[rj, rj, B 8 (x, y)]/5B e (x, y)\jj =v=0 = and is given 
by 

B e (x -y)= g 2 D(x - y)tr lC [A e G {x - y)]. (11) 

where Go stands for G[-Bq] and the trace is to be taken in Dirac and color space, whereas the 
flavor trace has been separated out. 

We will calculate the induced QCD vacuum condensates from the saddle-point expansion, 
that is, we will work at the mean field level. This is consistent with the large N c limit in the 
quark fields for a given model gluon two-point function. Under the mean field approximation, 
the dressed quark propagator G{x — y) is substituted by Go(x — y) which has the decomposition 

G \p) = i~f P + Z(p) = i 1 -pA(p 2 ) + B{p 2 ). (12) 
The X(j>) stands for the dressing self-energy of quarks and is defined as 

E(p) = A 9 Bl{p) = i 7 • p[A(p 2 ) - 1] + B(p 2 ), (13) 

where the self-energy functions A(p 2 ) and B(p 2 ) are determined by the rainbow Dyson-Schwinger 
equations (DSEs) 

t^ 2 > - * 2 - I / Wf s ' Dip - <14> 

B(P>) = % I (15) 



3 J (2^) 4y ^ HJ q 2 A 2 (q 2 ) + B 2 (q 2 ) 

Because the form of the gluon propagator g 2 D(s) in the infrared region is unknown, one often 
uses various model forms |15 | I16 | IT7] as input parameters in the previous studies of the Rainbow 
DSE. 

Within the DSE formalism, there are two qualitatively distinct solutions in Eg Q14|) and (|15(1 . 
The "Nambu-Goldstone" solution characterized by B(p 2 ) ^ describes a phase : (1) chiral 
symmetry is dynamically broken for it provides a momentum dependent constituent quark 
mass M(p 2 ) = B(p 2 )/A(p 2 ); and (2) the dressed quarks are confined for the dressed quark 
propagator does not have a Lehmann representation [TB] . The alternative "Wigner" solution 
characterized by B(p 2 ) = describes a phase with neither dynamical chiral symmetry breaking 
nor confinement. In this letter, we refer to the quark propagator in terms of the trivial solution 
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of the gap equation as the "perturbative" quark propagator Gq CT \p 2 ) with only the vector part 

[^ 2 >-*S/lF^- 9) ra>- (16> 

This "perturbative" quark propagator can be seen as the expectation value of the operator 
T[qi(x)qj(y)] over the perturbative vacuum \P) at the mean field level in the framework of 
GCM. 

From the GCM generating functional, it is now straightforward to calculate the vacuum 
expectation value(VEV) of any quark operator of the form 

O n ^ (^Ag lfe )% 2 Ag 2% ) • • • %, A£UJ, (17) 

in the mean field vacuum. Here the A® stands for an operator in Dirac, flavor, and color space. 
Take the appropriate number of derivatives with respect to external source terms rji and fjj of 
Eq. and set rji = fjj = [22] , we can get 

{o n ) = {-it D-mgl • • • A 51( G °U (1 ) • • • (GoW. (is) 

p 

where p stands for a permutation of the n indices. Once the dressing quark propagator Go(q 2 ) 
(We ignore the subscript below)is determined, one can calculate the two quark condensate 
(qq) , the four quark condensate^A^ggA^g) , etc. in the mean field level. Since the functional 
integration over the gluon field is quadratic in the framework of GCM, one can perform 
the integration over gluon field analytically. Using the same shorthand notation for the typical 
Gaussian integrations as in Ref. (Hj, we have 

J VAe~^ AD ~ lA+jA = e^ Dj 

j VAAe-^ AD ~ lA+ i A = (jD) e y Dj (19) 

J VAA 2 e- l 2 AD ~ lA+ i A = [D + {jD) 2 ]e^ 

where D is the dressing gluon propagator and is the quark color current. Because the 
gluon vacuum average can be replaced by a quark color current q^y^-^-q together with the 
gluon two-point function D, one can perform the integration over the quark operators in the 
mean field vacuum as described above. In this way, one can in principle obtain the vacuum 
expectation of value for any gluonic fields. This technique provides an feasible way to evaluate 
the expectation value of the operators with low-dimensional gluon fields such as the mixed 
quark-gluon condensate g^qG^a^q) in GCM (Note that for the VEV of operators with high 
powers of gluonic fields, this procedure will get rather complex). 
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3. Mixed Tensor Susceptibility 

With above preparation, the mixed tensor susceptibility can be calculated in the mean field 
level within this DSE model. The induced tensor susceptibilities %i K an d C are defined through 

(V\q^uq\V)z = 9 q xZ^(qq) (20) 
(V\qg c ^G« v q\V) z = g q nZ^{qq) (21) 
(V\qg c l5G^q\V}z = -ig q (Z^(qq), (22) 

where Z^ u stands for the external field, (V\ • • • \ V)z denotes the VEV over the QCD vacuum 
at the presence of the external filed and G^ = We^apG^ 0, . The nonzero VEVs of the 
operators above are due to the breakdown of Lorentz invariance in the presence of external 
constant field Zn V . From the QCD partition function for quarks in Euclidean space in the 
presence of the external field , the formulae for evaluating these susceptibilities take the form 

X(qq) = \J d 4 x(y|T[g(x)a^g(x),^ g ]|y) = ±IL x (0) (23) 

«(3?> = \ j d A x{V\T[q{x)g c ^G%q{x),q<j, lu q]\V) = ^U K (0), (24) 

according to Ref. [HI El- Due to the fact that the vacuum susceptibilities reflect the nonper- 
turbative structure of the QCD vacuum, n x (0) and II K (0) on the right hand side of above 
Eas. (|2*3*|) and (|2*1)) should be subtracted by the corresponding perturbative contribution terms. 
Within the DSE formalism, the perturbative contribution to n x (0) and II K (0) can be evaluated 
by the trivial quark propagator, namely the "perturbative" quark propagator in terms of the 
trivial Wigner solution to the dressed quark gap equations (|14(l and (|15|) . This is a reasonable 
procedure because the Wigner solution of the dressed quark propagator describes a phase with 
neither DCSB nor confinement and the difference between the Nambu solution and the Wigner 
solution vanishes at short distance according to numerical studies. In fact, the Wigner solution 
can play the role of the perturbative dressed quark propagator has been used extensively in the 
study of thermal property of QCD within the DSE formalism^?], where the bag constant was 
defined as the difference of pressure between the true QCD vacuum and the perturbative QCD 
vacuum, which were evaluated by the Nambu-Goldstone solution and the Wigner solution to 
the quark propagator, respectively 18 . 
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Therefore, we rewrite the Eqs. p3|) and (|2i]> as 

1 



X(qq) = § I d 4 x(V\T[q(x)a^q(x),qa^q}\V) N 



~ J d A x{P\T[q{x)o llv q{x),qa t _ lu q}\P) 



^7(0), (25) 



«<??) = ^^(yir^^yG^^x),^^]^)^ 

- l - J d A x{P\T[q{x)g c ^-G%q{x),qa, u q]\P 



= ^7(0). (26) 

By substituting the "perturbative" quark propagator G per (p 2 ) to Eq.(|18|). the determination of 
the expectation value of the T-product operators in terms of quark fields over the perturbative 
vacuum state \P) can be performed self-consistently within the GCM formalism. It should be 
noted that the evaluation of x{qq) i n Ref. 11 is consistent with Ea. ((2*5|) because in this special 
case the subtraction terms to Eq. (|23)) has zero contribution to x(qq) due to B'(p 2 ) = . 

Using Eq. (|T!I|) , the expression for VEV of above T-product operator including gluonic fields 
is converted to the VEV for the product with the form of (|17|) in terms of six quark fields and 
eight quark fields. According to Eq. (|18|) . we have 



i J d A x(V\T[q(x)g c ^G a ^),q^uq}\V) N = 

-|i J dx 4 J dz 4 g 2 [d;D(z-x)] x 

ltr D [G(x - z) lv G{z - Q)a lw G(0 - x)] 

+tr D [G(x - 0)<vG(0 - z) lv G(z - x)] X 

-2i J dx A J dz\ J dzjg 2 D( Zl -x)g 2 D(z 2 -x) x 

itr D [G{x - zxYfpGizx - z 2 ) lu G(z 2 - 0)^(2(0 - x)] 
+tr D [G(x - zi)7 M G(*i - 0)a IMJ G(0 - z 2 ) lv G{z 2 - x)] 
+tr D [G(x - 0)a flu G(0 - z^Gfa - z 2 ) lv G{z 2 -x)}\. (27) 

Substituting G per (x — y) for G(x — y), the similar expression for the VEV of same operator 
over the perturbative vacuum can be obtained. After Fourier transformation, we find the first 
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part of right hand side of (|27|) is zero and the final result for the mixed tensor susceptibility in 
the momentum space takes the form 

«<3ff) = T ^Jdss[^] 2 [fB\s) + ^-sA(s)(2-A(s))] 

^ J dsdt J dxsty/l - x 2 g 2 D(s,t,Vsix)Z- 2 (s)Z- 1 (t) 



+ 



32 

:B(s)^B{s)B 2 (t) + B(t)A(s)A(t)Vsix - [A{s) - 1 

; [2A{s)B{t)Vsix - A{t)B{s)] 
3 

32vr 5 



J dsdt J dxst^/l - x 2 g 2 D(s, t, Vsix)Z'- 2 (s)Z'- 1 (t) 
xA' 2 (s)A'(t) [A'(t) - 1] [st - Astx 2 ] 

^-r f dsdt [ dxst^/T^x^g 2 D(s,t,Vsix)Z- 2 (s)Z- 1 (t) 

327T & J J_i 

xA 2 (s)A(t) [A(t) - 1] [st - Astx 2 ] , (28) 

where Z{s) = sA 2 (s) + B(s) and Z'(s) = sA' 2 (s). 

It should be noted that to get this expression the Dyson- Schwinger equation (fTTj) has been 
used again. The UV divergence of Eq. (|27f) can be illustrated by a simple analytical confining 
model g 2 D(p — q) = -^(2ir) 4 r] 2 d 4 (p — a) |19j. In this model, the expression for (|2*7|) takes a 
relative simple form 

1 f d 4 x(V\T[q(x)g c ^-G%q(x),q^ u q]\V) N = 



6 



h I dss[^] 2 [^B\s) + 2 lsA{s){2-A{s))] 
dssZ- 3 (s) [B 4 (s) + sB 2 (s)A(s)} 



167T 

36r/ 2 
16vr 2 

„2 



+ S / dss3z ~ 3 ( s )i A ( s ) ~ l]A*)- (29) 



The Nambu-Goldstone solution for this model is 



2 v 2 < rj - 

Mp 2 ) = < , „ , rn?, , 4 ' (30) 



2 (1 + W1 + ^-), otherwise, 



B{p 2 ) = { v ' " ' 4 ' (31) 

0, otherwise. 
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The alternative Wigner solution takes the form 

BV) = 0, A'(p 2 ) = ^l + Jl + ^ r ). (32) 

Due to A(s) — 1 ~ c/s for s — > oo according to . the last term of right hand side of Eq. 
(|29|) is logarithmic divergent. In addition, replacing A(s) with A'(s) to (|29|). there still exits 
logarithmic divergence due to -A'(s) having the same behavior as -A(s) in the large energy region. 
Because the vector part A(s) as well as the scalar part B(s) both reflect the nonperturbative 
information in the low energy region, it is more reasonable to subtract the corresponding 
perturbative part rather than to simply ignore this divergent term in (|29j) . 

With the effective subtraction of the perturbative contribution, there is no UV divergence in 
above integrations. Actually, the subtraction procedure guarantee K(qq) the role of the order 
parameter for QCD chiral phase transition because it becomes zero when QCD undergoes a 
phase transition from the Nambu-Goldstone phase to Wigner phase (That means A(p 2 ) — ► 
A'{p 2 ) and B{p 2 ) -»0). 

4. Results and Discussions 

The determination of the mixed tensor susceptibility is based on the same effective gluon 
propagator models g 2 D(s) which had been used in Ref. ^?J[2U|. In general, the quark-quark 
interaction g 2 D(s) has the form 

9 2 D(S) = (33) 

S 

where s = p 2 . Two popular quark-quark interaction models with two parameters for a(s) are 
investigated here: 

ai(s) = 3vrsx 2 ^2 (34) 

a 2 (-s) = d-rrs X - , (35) 
s z + A 

where d = fX. The two low-momentum parameters, the stength parameter x an d the range 
parapmeter A , are varied with the pion decay constant fixed at 87 MeV which is more ap- 
propriate in the chiral limit rather than the pion's mass-shell value of 93 MeV. Noted that the 
above quark-quark interactions dominate for small s and simulate the infrared enhancement 
and confinement. Because the effective quark-quark interactions (|34I35|) have a finite range in 
momentum space, the momentum integral for the calculation of the quark condensate 

3 f°° r B (s) , x 

m = -\^ 2 h dss sA 2 ( s ) + B 2 ( s y (36) 

is finite. According to the obtained values of the chiral low energy coefficients Li following 
Ref. [JO] based on both ansatz (|34j) and (|35|) are compatible with the phenomenological values. 
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The model ansatz (|35|) has been successfully used to investigate the space structure of the 
non-local quark condensate (q(x)q(0)) in Ref.[2l. within GCM formalism. It should be stressed 
in this context that our interactions are not renormalizable due to using the bare quark-gluon 
vertex within the rainbow DSE formalism. Therefore, the scale at which a condensate is defined 
in our calculation is a typical hadronic scale, which is implicitly determined by the model quark- 
quark interaction and the solutions of the DSEs for the dressed quark propagator. The similar 
case is the determination of the vacuum condensate in the instanton liquid model [23] where 
the scale is set by the inverse instanton size. 

To check the sensitivity of the mixed tensor susceptibility on the forms of quark-quark in- 
teractions, the above models with different sets of parameters x an d A are investigated below, 
where the results for the quark condensate are also given. 

Table 1. The value of n{qq) for model 1 with two sets 



of different parameters. 



A(GeV 2 ) 


x(GeV) -(qq)-s(MeV) 


K{qq)(GeV 4 ) 


re(GeV) 


0.2 


1.55 213 


2.5* 10" 3 


-0.26 


0.02 


1.39 170 


2.9 * 10~ 3 


-0.59 


0.002 


1.30 149 


1.4* 10" 3 


-0.41 


Table 2. 


The value of n{qq) for model 2 with four sets 




of different parameters. 




A(GeV 4 ) 


x(GeV) -(qq)i(MeV) 


K(qq)(GeV 4 ) 


re(GeV) 


1 * 10" 1 


1.77 290 


4.7* 10" 3 


-0.19 


1 * 10~ 2 


1.33 250 


3.6* 10" 3 


-0.23 


1 * 10~ 4 


0.95 217 


3.0* 10" 3 


-0.29 


1 * 10~ 6 


0.77 204 


3.1 * 10~ 3 


-0.36 



In Table 1 we display the values for (qq) and n(qq) based on model 1 with three sets of 
parameters and in Table 2 the same quantities with four sets of parameters based on model 2. 
In both cases, the obtained values for {qq) are compatible with the standard phenomenological 
value in SVZ sum rules, whereas the mixed tensor susceptibility K,(qq) is much small. The 
previous estimation of n(qq) obtained In Ref. |Hj is 0.10 GeV 4 . Actually, the value of the quark 
condensate tensor susceptibility x(qq) obtained within GCM formalism is also very small 
compared with the estimation based on SVZ sum rules . In fact, different versions of SVZ 
sum rules have given very different values for x(qq) m previous studies (HI El El- Therefore, it 
shows that the induced vacuum condensates have very little impact on the determination of 



10 



the nucleon tensor charge from the theoretical formalism of DSEs. 

In summary, we have investigated the mixed tensor susceptibility at the mean field level in 
the framework of GCM/DSE formalism. In the calculations, the vacuum matrix elements for 
the operator in terms of quark and gluonic fields can be obtained by substituting the gluonic 
fields with the quark color current operator and the model gluon propagator which describes 
the effective quark-quark interaction within the GCM formalism. To subtract the perturbative 
contribution to the expression for the mixed tensor susceptibility, the Wigner solution to the 
quark gap equation was used self-consistently in this formalism. Using different quark-quark 
interaction models, we find that the mixed tensor susceptibility as well as the quark condensate 
tensor are both very small susceptibility does within DSE formalism. Therefore, we get the 
conclusion that the induced vacuum condensates have little effect on the determination of the 
nucleon tensor charge. Finally, we want to tress that above approach can also be used to 
investigate the other mixed susceptibility of the QCD vacuum. 
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